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NEW  METHODS  OF  MEASURING  NORMAL  ACOUSTIC  IMPEDANCE 

James  L.  Wayman 

Department  of  Mathematics 
U.S.  Naval  Postgraduate  School 
Monterey,  Calif.   93940 

INTRODUCTION 

About  three  years  ago,  owing  to  our  work  in  computer  simulation 
of  sound  fields  in  enclosures,  we  developed  the  need  for  closely 
spaced,  broad-band  measurements  of  normal  acoustic  impedance  for  a 
variety  of  architectural  materials.  Through  text  books  and  student 
laboratories  I  was  well  aware  of  the  Standing-Wave-Ratio  technique 
and  the  tedium  of  its  application.  Given  our  problem  and  the  extreme 
numbers  of  data  points  we  needed  to  collect,  a  technique  orders  of 
magnitude  faster  was  clearly  required.  An  answer  was  found  in  a 
1977  paper  by  A.F.  Seybert  and  D.F.  Ross  (Ref.5).  During  the  course 
of  our  research  with  the  Seybert/Ross  technique,  another  related 
method  was  published  in  the  Journal  by  J.Y.  Chung  and  D.F.  Blaser 
(Ref.2).  We  found  the  Chung/Blaser  technique  to  be  as  easy  to 
implement  and  as  accurate,  but  to  compute  in  a  much  shorter  time 
than  the  Seybert/Ross  method.  The  purpose  of  this  paper  will  be 
to  discuss  and  compare  the  implementation  of  these  two  related  methods. 

THE  HISTORY  OF  DUAL,  FIXED-MICROPHONE  MEASURING  TECHNIQUES 

It  surprised  me  to  learn  that  dual,  fixed-microphone  techniques 
for  measuring  the  normal  acoustic  impedance  of  materials  are  not 
new,  but  date  back  as  far  perhaps  as  1932.  The  first  mention  in  the 
literature  was  a  1941  paper  by  Clapp  and  Firestone  on  acoustic 
wattmeters  (Ref.3).  They  recognized  that  their  device  could  be  used 
to  measure  the  absorption  coefficients  of  materials.  Then  at  the  28th 
meeting  of  the  Acoustical  Society  in  1943  R.H.  Bolt  and  A. A.  Petrauskas 
delivered  a  paper  titled  "An  Acoustic  Meter  for  Rapid  Field  Measurements 
(Ref.l).  Some  years  later,  Dr.  Ted  Schultz  wrote  his  Ph.D.  thesis 
and  published  a  paper  on  the  acoustic  wattmeter  and  discussed  its 
application  to  the  measurement  of  absorption  coefficients  (Ref.6). 
These  methods  all  used  analog  circuitry  and  sinusodial  sound  sources. 
Perhaps  one  of  the  reasons  that  use  of  the  acoustic  wattmeter  for 
impedance  measurements  has  beneclipsed  by  the  Standing-Wave-Ratio 
apparatus  has  been  the  complexity  of  the  wattmeters  circuitry. 

Recently  S.J.  Elliott  has  published  a  description  of  a  low-cost, 


simple  system  of  this  type  (Ref.4  )  so  perhaps  the  use  of  dual,  microphone 
analog  techniques  for  impedance  measurement  will  gain  popularity  in  the 
future. 

DIGITAL  TECHNIQUES 

Both  the  Seybert/Ross  and  the  Chung/Blaser  techniques  to  be 
discussed  today  are  digital  methods  and  consequently  require  digital 
signal  processing  equipment.  Figure  1  is  a  schematic  diagram  of  the 
necessary  equipment.  The  core  pieces  of  equipment  are  a  tube,  a 
desktop  computer,  and  a  signal  analyzer.  A  speaker  one  end  of  the 
tube  is  excited  by  white  noise  and,  with  an  unknown  impedance  at 
the  other  end  of  the  tube,  the  auto-  and  cross-spectra  of  the  the 
signal  at  the  two  microphone  locations  is  measured.  From  these 
values  and  the  knowledge  of  the  microphone  spacings,  the  acoustic 
impedance  of  the  material  at  the  end  of  the  tube  can  be  calculated. 
Unlike  the  Standing-Wave-Ratio  method,  these  techniques  do  not  require 
that  the  sample  material  be  placed  inside  the  tube.  Rather,  the  tube  is 
placed  against  the  sample,  thus  allowing  in  situ   measurements. 

THE  SEYBERT/ROSS  METHOD 

Let's  look  first  at  the  principles  of  the  Seybert/Ross  method 
because  mathematically,  this  is  the  simplier  of  the  two.   The 
frequency  dependent  absorption  coefficient  a(f)  of  a  material  is  unity  minus 
the  ratio  of  reflected  sound  power  to  the  incident  sound  power.  This 
can  be  given  by  c  f£\ 


a(f,)=  l  '  -5^W 


where  Sry,(f)  and  S,.j  (f)  are  the  frequency  dependent  auto-power  spectral 
values  of  the 'reflected  and  incident  waves.  Similarly,  the  phase 
change  upon  reflection  can  be  given  by 


<j>(f)=  tan"1 


where  Qir  and  Cir  are  the  imaginary  and  real  parts  of  the  cross-power 
spectrum  of  the  incident  to  reflected  waves. 

Note  that  we  cannot  measure  S.^f),  Srr(f),C1r(f)  and  Qi  (f) 
directly,  but  we  can  measure  Su(f),  S22(f),  C12(f)  and  Q12(f),  the 
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auto-  and  cross^power  spectra  of  the  signals  at  microphones  1  and  2.  The 
signals  at  microphones  1  and  2  contain  both  incident  and  reflected  waves. 
Can  S,,(f)  etc.  be  decomposed  to  yield  S..(f)  etc.? 

The  answer  is  "yes"  and  quite  readily  in  fact.  The  mathematics  is 
marvel ously  simple,  requiring  only  the  definition  of  a  finite  Fourier 
transform,  Euler's  relations  and  elementary  complex  algebra.  I'll  spare 
you  the  details  (they're  contained  in  Ref.5  ),  but  with  an  additional  ten 
minutes,  I'm  sure  I  could  convince  you  that 

Sn(f)  =  S^f)  +  Srr(f)  +  2  cos  2kX;L  C-r(f)  +  2  sin  2kX;L  Qir 

where  k  is  the  wave  number  and  x.  is  the  distance  between  microphone 
1  and  the  sample  end  of  the  tube.   Note  that  S,.(f)  is  a  linear  function 
of  S^f),  Srr(f),  Cir(f)  and  Q1p(f),  S22(f),  C12(f)  and  Q12(f)  are  as 
well.  Thus  we  can  write 


sn(f) 
s22(f) 
c12(f) 
Q12(f) 


[aJ 


S^f) 

srr(f) 
c  (f) 
Q.  (f) 


where  £aJ  is  the  coefficient  matrix.  The  good  news  is  that  [Ajis,  in 
general,  invertable.  The  bad  news  is  that  the  coefficients  are  functions 
of  the  wave  number  k  ,  making  [a] frequency  dependent.  Consequently,  [aJ 
must  be  inverted  at  each  frequency  of  interest  and  multiplied  by  the 
vector  on  the  left-hand  side  to  yield  the  auto-  and  cross-power  spectral 
values  for  the  incident  and  reflected  waves.  From  these  values,  then,  the 
absorption  coefficients  and  the  phase  change  can  be  computed  as  previously 
mentioned. 

If  our  orginal  goal  had  been  to  find  the  complex  acoustic  impedance 
rather  that  the  absorption  coefficient  and  phase  change,  we  would  have 

started  with  the  relation 

Cn..(f)  +  i(L.(f) 
Z(f)  = 


M 


£u 


S  (f) 
uuv  ' 


where  Z  is  the  complex  acoustic  impedance,  p  is  the  acoustic  pressure,  u 
is  the  particle  velocity,  and  the  other  symbols  are  defined  as  before. 
Using  standard  acoustic  relationships,  we  could  find  Z(f)  in  terms  of 
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SlM(f),  Srr(f),  Cir(f)  and  Qir(f). 

THE  CONSTRUCTION  OF  THE  TUBE 

Before  we  consider  the  Chung/Blaser  technique,  let's  consider  the 
details  of  the  impedance  tube.  We  wanted  to  make  our  first  tube  as 
quickly  and  as  cheaply  as  possible.  Essentially  we  took  a  3"  round 
speaker  and  bolted  it  to  one  end  of  a  3"  diameter,  12"  long  piece  of 
PVC  pipe.  Two  holes  were  drilled  into  the  tube  wall  and  some  Tygon 
tubing  inserted  to  hold  the  7/8"  diameter  microphones  in  place. 

The  HP-5420  signal  analyzer  has  its  own  low-pass  filtered  white  noise 
source.  The  noise  must  be  filtered  to  remove  frequencies  above  those 
we  are  measuring,  and  thereby  prevent  aliasing  phenomena.   The  noise 
was  amplified  and  sent  to  the  speaker.  The  microphone  signals  were 
preamplified  and  sent  to  the  HP-5420  where  auto-  and  cross-spectral  data 
was  averaged  over  30  samples  and  resolved  into  256  frequency  bins. 
This  data  was  dumped  into  an  HP-85  controller  where  it  was  multiplied  by 
the  inverted  coefficient  matrix  calculated  at  each  frequency  as  discussed 
above.  Results  were  plotted  from  the  HP-85.  The  entire  process  of  data 
collection,  calculation  and  plotting  took  roughly  18  minutes.  We  were 
able  to  obtain  data  from  about  200  Hz  to  2500  Hz.  The  lower  limit  seems 
regulated  by  our  ability  to  measure  phase  differences  between  the  microphone 
signals  for  wery   large  wavelength.    The  upper  limit  is  controlled  either 
by  the  microphone  spacing,  which  must  be  less  than  h  wavelength,  or  by 
the  first  "sloshing"  mode  of  the  tube,  above  which  our  assumptions  of 
plane  waves  travelling  down  the  tube  and  reflecting  break  down. 

MICROPHONE  CALIBRATION 

We  must  mention  here  the  phase  and  amplitude  calibration  of  the 
two  microphone/preamplifier  sections.  There  are  basically  two  techniques 
available  (and  perhaps  some  creative  combinations  of  the  two).  We  could 
stop  the  data  collection  half  way  through  and  switch  microphone  systems,  or 
we  could  place  the  microphones  in  identical  sound  fields  before  we  start  and 
record  any  frequency-dependent  phase  and  amplitude  differences.  Test 
data  could  then  be  preconditioned  in  the  controller  by  this  record  before 
multiplication  by  the  inverted  matrices.  Both  methods  are  discussed  in 
the  literature.  We  found  that  the  latter  method  was  superior  as  it  allowed 
for  imbalances  in  the  A/D  converter  sections  of  the  signal  processor,  whereas 
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the  microphone  switching  method  did  not.  The  identical  sound  fields  were 
obtained  by  placing  the  microphones  in  a  plate  mounted  at  the  sample  end 
of  the  tube,  assuming  the  sound  field  to  be  radially  symmetric  at  frequencies 
below  sloshing. 

VALIDATION  OF  TUBE  PERFORMANCE 

The  performance  of  the  tube  was  then  validated  against  a  known  system. 
As  suggested  in  the  Seybert/Ross  paper,  we  used  a  piece  of  pipe, the  same 
diameter  as  the  impedance  tube,  with  a  reflective  cap  to  provide  a 
purely  reactive  impedance.  We  would  expect  a  absorption  coefficient  of 

0  and  a  phase  change  equal  to   «i 

—  .~2* 

Good  results  were  obtained  out  to  2500  Hz. 

THE  CHUNG/BLASER  METHOD 

We  wanted  to  improve  the  system  in  two  ways:  we  wanted  a  higher 
frequency  range  and  we  wanted  to  cut  down  the  18  minute  running  time. 
We  changed  the  system  controller  to  a  faster  HP-87,  built  a  fancy 
aluminum  tube  of  IV  diameter  and  used  V  diameter  microphones  with  a 
closer  spacing.  Simultaneously,  we  discovered  the  1981  Chung/Blaser 
paper.  The  system  described  by  Chung  and  Blaser  is  the  same,  but  the 
mathematical  formulation  is  different. 

Chung  and  Blaser  show  that  the  complex  reflection  coefficient 

can  be  expressed  in  terms  of  acoustic  transfer  functions  between  the 

two  microphones  as 

H12(f)  -  H^f) 


R(f)  = 


Hr(f)  -  H12(f) 


where  H19(f)  is  the  acoustic  transfer  function  between  microphone 

C  +iO 
1  and  2  and  is  expressible  by    12  v12 

and  H.  is  the  transfer  function  of  the  incident  wave  between  the  two 
microphone  locations  and  H  (f)  is  the  transfer  function  of  the  reflected 
wave  between  the  microphone  locations.  The  transfer  functions  of  the 

incident  and  reflected  waves  between  the  microphones  are  simply  the 

•1  ks      i  ks 
phase  delays  associated  with  the  microphone  spacing  e    and  e 
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where  s  is  the  microphone  spacing.  Be  careful  to  note  that  the  microphones 
in  the  Chung/Blaser  paper  are  numbered  oppositely  from  the  Seybert/Ross 
work. 

The  acoustic  impedance,  absorption  coefficient,  and  phase  change 
can  all  be  readily  computed  from  the  complex  reflection  coefficient. 
The  above  equations  showed  that  the  complex  reflection  coefficient  could 
be  calculated  from  the  cross-power  spectra  of  the  two  microphones  and 
the  auto-power  spectra  of  microphone  1  alone,  without  matrix  inversion. 
Therefore,  this  method  should  save  us  time  in  both  data  transfer  from 
the  signal  processor  and  in  calculation  time.  Using  this  new  system, 
in  fact,  running  time  was  cut  to  about  12  minutes. 

RESULTS  USING  THE  NEW  SYSTEM 

The  system  was  again  validated  against  a  capped  tube.  Results 
showing  experimental  and  theoretical  results  are  attached.  Good  agreement 
is  obtained  out  to  about  4000  Hz.  Also  attached  are  results  obtained  for 
several  other  architectural  materials. 

COMPARISON  OF  SEYBERT/ROSS  AND  CHUNG/BLASER  RESULTS 

We  wondered  if  there  were  any  computational  differences  affecting 
accuracy  between  the  two  methods.  Both  methods  were  put  to  work  on 
the  same  auto-  and  cross-spectra  data  sets.  The  results,  which  are  attached, 
appear  identical,  although  the  Seybert/Ross  system  required  240  seconds 
of  computing  time  versus  100  seconds  for  the  Chung/Blaser  method. 

CONCLUSIONS 

We  have  been  very  pleased  with  our  two  systems,  but  now  use  the 
Chung/Blaser  method  exclusively  because  of  its  faster  processing  time. 
Output  of  the  method  can  be  stored  in  magnetic  form  for  other  uses,  which 
is  extremely  helpful.  We  are  building  an  even  smaller  tube,  using  smaller 
microphones  with  still  closer  spacing,  hoping  to  push  our  measurement 
capability  to  6400  Hz  and  beyond. 
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Figure    1.    Experimental   Set-up   Simplified   Block   Diagram 
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